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, We introduce an generalized action functional describing the equations of motion and the vari- 

' ational equations for any Lagrangian system. Using this novel scheme we are able to generalize 

I Noether's theorem in such a way that to any n-parameter continuous symmetry group of the La- 

fSJ ■ grangian there exist 1) the usual n constants of motion and 2) n extra constants valid in the vari- 

ational equations. The new constants are related to the infinitesimal generators of the symmetry 
, transformation by relations similar to the ones that stem from the 'nonextended' Noether theorem. 

m 

The equations of motion of many physically relevant systems can be obtained within a Lagrangian framework. 
The systems are described using a Lagrangian density through the Euler-Lagrange equations. To these equations 
of motion we can associate another set of equations obtained by linearizing the Euler-Lagrange equations around 
' one of its particular solutions. Such linearized equations are called the variational equations [1], are important in 
\^ ' studies of stability [2-4] and for defining quantities as the Lyapunov exponents or some forms of entropy in dynamical 
system theory [2,5]. Also they are the equations of geodesic deviation in general relativity and other metric theories 
' of gravitation [7-11]. They are in general useful for describing perturbations from the original dynamics and have 
, been variously used in analysing non-linear evolution equations [12], in studies of stability in galactic dynamics [13], 
in performing the Painleve test [6], and in geometric control theory [14,15]. 

It has been shown that if the field equations can be casted in Hamiltonian form and if any solution of the variational 
equations can be expressed as the Poisson bracket of a field variable, i/ia, with any other dynamical variable, then this 
dynamical variable is necessarily a constant of motion [12]. The converse of this result has been proved earlier in [16]. 
The variational equations and its symmetry properties are also important for studying solitonic solutions to certain 
non-linear wave equations [4] . It is also known that if the original Lagrangian is invariant under a translation in any 
space-time direction and ip is one of its solutions, then the directional derivative of ip along such direction is then a 
solution of the variational equations [16]. 

Thus, given the relevance of the variational equations in many aspects of physics, and of the relationship [12,16] 
between constants of motion and properties of the solutions of the original equations, in this work we want to pursue 
. ^ the relation between symmetries and constants of motion in the first-order dynamics. To accomplish this task, a 
■ Lagrangian description of the variational equations is very convenient. Though this description does not appear easy 
" " ' to attain — given the use of particular (and explicitly space-time dependent) solutions of the original equations for 
deriving the variational equations — a complete Lagrangian characterization of the variational equations of Lagrangian 
systems has been advanced recently [17], see also [8,11] for previous related results. The idea is the following, given 
the Lagrangian density associated with a system, L, a new density needs to be defined [14,17] as 

dL ^ dL ^ 

7 = + (1) 



o 



X 



where the and the correspond, respectively, to deviations from the original field variables and to their space-time 
derivatives. These deviations are assumed to connect two nearby solutions of the original Euler-Lagrange equations. 
In (1), as in all of the paper, the summation convention is implied, we also use ip^f^ = dijj/dx^ , Latin indices correspond 
to internal variables, and Greek indices to space-time coordinates: a;**, = 0, . . . , 3. From a strictly mathematical 
viewpoint the density 7 can be interpreted as the prolongation of L [14]. 
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Using the density 7 we can further define an alternative action functional as 



E = 



where Q, is an appropriate space-time region. This action is fundamental for our formulation. 
The extremalization of (2) gives [17] both the Eulcr-Lagrange equations, 



dL 



dL 



0, 



(2) 



(3) 



and the associated variational equations, 



0. 



(4) 



In the form given above, the variational equations do not have to be regarded as explicitly dependent on the space- 
time coordinates unless the original Lagrangian density L is so from the start. This property, shared by the function 
7, is analogous to what occurs with the Lagrangian density L itself, which is not considered space-time dependent 
despite being an explicit function of the fields which, in principle, also depend on x^. The form (4) of the 

variational equations and their origin from (2) is closely connected with properties of the Jacobi equation of interest 
for the study of geodesies in Riemannian manifolds and gives a way for deriving the associated curvature tensor from 
a variational formulation [2,18]. It should be clear that 7 plays the role of a new form of the density associated with 
the system. From the existence of the action (2) it can be easily proved that 7 is invariant under arbitrary point 
transformations on its configuration manifold comprised by the (j)'^, the e", and their space-time derivatives [14,17]. 
The function 7 thus describes both the system and its variational equations. 

To establish the connection between symmetries and conservation laws in the variational equations, let us consider 
a prolonged Lagrangian density 7 and assume it invariant (or quasi-invariant, but in what follows we use invariant 
for short) under the following continuous group of transformations 



(5) 
(6) 
(7) 



where the w'^ are the r parameters of the group, defined in such a way that the transformation reduces to the identity 
when all the w" vanish. The infinitesimal generators of the transformation group (5) are thus 




(8) 



If the set of transformations (5 7) [or (8)] is a symmetry of 7, it could stem directly from a symmetry of the original 
Lagrangian L and, in this case, we could have = F^. But it is also possible that this will not be the case, for 
the symmetries of 7 are expected to be larger than those of L — as can be simply ascertained by thinking about the 
possible transformations that leave (3) and (4) unchanged. Our examples below show simple instances of such case. 

That the transformation group (5) be a symmetry of 7, means that the action S (Eq. 2) remains invariant under 
the transformations generated by the group. A necessary and sufficient condition for this to be true is that 



7 



D : 



(9) 



where D = det (dx^ / dx") is the Jacobian determinant of the transformation. Eq. (9) is the fundamental equation 
allowing us to explicitly write the relationship between the generators of the continuous symmetries of 7 and the 
conserved quantities in the variational equations [19]. 

To see how the extension can be accomplished, differentiate partially this equality (9) with respect to w^, noting 
that the right hand side is independent of such parameters, and evaluate in = after the differentiation. After 
some manipulations this yields 
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+ 
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( dl ^ 







(c-ec) 

'Til - = 0. 



Finally, use the equations of motion [Eqs. (3) and (4)] to obtain 



^=0 



where the divergenceless tensor rj* is defined by 



(10) 



(11) 



(12) 



This is the conserved tensor associated with the variational equations. It is to be noted that this result can be 
regarded as the Noether's theorem for the variational equations that can be checked following the steps for the 
demonstration of Noether's theorem that can be found in many places, for example in [19]. The conservation of r 
(12) was to be expected since 7 is indeed a Lagrangian function on its own. The conserved quantities are simply 
related to the infinitesimal generators (8) of the symmetry group of the variational equations, as (12) explicitly 
shows. Notice that our results also show that in any Lagrangian system invariant under a n-parameter continuous 
group of transformations there exist 2n conserved quantities, the usual n Noether conserved quantities plus the n 
additional ones associated with the variational equations. Of course that not very one of the n new quantities has to 
be independent of the n previously known ones. 

To give examples of the usefulness of this result, we next derive (in admittedly simplistic instances) conserved 
quantities starting from well-known Lagrangian symmetries. 
Example 1. Let us consider the Lagrangian density 



2^ ^ '^a,ju^6,i^5 

where i?"^^'' = i?/"""^ is the Riemann tensor. 

This Lagrangian (13) describes Einstein equations in a vacuum. The density 7 associated with (13) is 



(13) 



1 



7 



(14) 



in this case 7 describes linearized gravitation in a vacuum. Taking into consideration the properties of R"^^"^, it is 
not difficult to realize that 7 is invariant under the one parameter group of transformations 



and 



(15) 



the infinitesimal generators are thus ^ = 0, = x'^ , and 77'* = 0. With the symmetries and the infinitesimal generators 
established, the extended Noether's theorem predicts that the quantity 
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{R''^''^)^^ecXt,,^Xa + R''^""'et„.Xa 



(16) 
(17) 



is divergenceless and hence generates a conserved quantity. This is a rather useful tensor that was employed (and had 
to be evaluated using direct calculations) in [20] to calculate the conserved quantities of Newman and Penrose. Notice 
also that this is a conserved quantity in any vacuum spacetime in General Relativity. We have thus managed to 
derive using our result (12) one of the results in [8]. Moreover, 7 can be also shown to be invariant under the one 
parameter group of transformations 



and 



(18) 



3 



with the infinitesimal generators ^ = 0, C** = 0, and t]^ = e^. In this case (12) leads directly to the divergence less 
tensor 



= e„ar6,., (19) 

which due to the symmetry properties of i?"^^"^, happens to identically vanish. Nevertheless, we have illustrated in 
a simple instance that the symmetries of 7 can be larger than those of L. It should bo noticed that the 'extended' 
symmetry [equations (18)] is the original symmetry of the Lagrangian just applied to the perturbation variables. We 
must point out also that the idea for this example is taken from [8] where it is used to relate symmetry transformations 
to the existence of divergence less quantities in a first-order approximation to Einstein theory. 

There are other potential uses of our main result (11) in general relativity and in other field theories. For example, 
for investigating perturbations to known metrics [21] and to contribute to the study of certain angular momentum 
ambiguities of recent interest [22] . The excellent review [23] discusses some other contemporary uses of perturbation 
methods in cosmological models in the context of general relativity. 

Example 2. Our result can be also profitably applied to particle mechanics by a simple reinterpretation of the symbols 
used in Eq. (12) and the replacement of the four spacetime parameters by the time t. To give a simple example, 
let us consider a Lagrangian function describing a system of particles in which one of the coordinates is ignorable. In 
this case, the original Lagrangian and the function 7 can be written as 

L = ^mabrq' - U{q), (20) 
fjjj 

7 = m„,g«e^ - — e^ (21) 

where niab is a symmetric 'mass' matrix formed by the second derivatives of the system's kinetic energy respect to the 

generalized velocities q", and U{q) is its potential energy function. For the purposes of the example we arc assuming 
that niab does not depend on the generalized coordinates g", but this is of no important for the argument — it just 
simplifies a little the form of the associated conserved quantity. The invariance under translations in the, let us say, 
q"^ direction (i e. the transformation is q^ = q"^ + w, — + w, the rest of the coordinates and velocities remain 
unchanged) can be associated with the infinitesimal generators ^ = and = t]^ = 1, and so using (12) we can 
obtain that 

de^ dq^ 



= P. + ^, (22) 



is a constant of motion in the variational equations of the system, where pA = dLjdqA is the momentum conjugate 
to the coordinate q^ . In fact, using the obvious invariance of (20) under the space translation q^ = q^ ^ w, and 
e'^ = e"* + w, which corresponds to the infinitesimal generators ("^ = 1, rj"^ = 1, and ^ = 0, we can directly prove that 
the term 



d'y 
dq 



(23) 



appearing in Eq. (22), is, as it should be, a constant of motion on its own. These results are direct and simple examples 
of conserved quantities which may be important (in more complex instances) in the analysis of classical perturbation 
methods with astronomical interest [24]. 

To summarize, we have established a form of Noether's theorem that encompasses the variational equations and 
have illustrated its use with direct examples. The conserved quantities encountered can be useful in perturbation 
theory as our formalism can be regarded as a starting point for studying perturbations in all kind of Lagrangian 
problems. Moreover, if we take into consideration the Lagrangian foundations of the path integral approach our 
formalism could be useful in approaching approximate schemes in quantum field theory [25]. We also consider it as 
possibly important for studying properties of solitonic solutions in nonlinear equations [4]. Besides, the formalism 
could have some bearings to modern mathematical developments in Lagrangian field theories [26]. 
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